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^ . Abstract 

a • 

^ ' A parametrization of density operators for bipartite quantum systems is proposed. It is based 

Q^l on the particular parametrization of the unitary group found recently by Jarlskog. It is expected 

that this parametrization will find interesting applications in the study of quantum properties of 
I , many partite systems. 

> ■ 

oo : 

(N 1 Introduction 
vn 

Density operators represent states of quantum systems. They are crucial to describe the dynamics 

^ \ of open quantum systems [H [2] . Recently, there is a considerable interest in the structure of the set 

(yQ of density operators due to the emerging field of quantum information theory [3]. It turns out that 

■ quantum entanglement may be used as basic resource in quantum information processing and com- 

J> . munication. The prominent examples are quantum cryptography, quantum teleportation, quantum 

'k^ [ error correction codes and quantum computation. It is, therefore, clear that the proper description of 

^ • the set of density operators is highly important. 



It turns out that the structure of the set of density operators is nontrivial and it is well understood 
only for 2-level systems [H [5] . This problem was studied by many authors from different perspectives 
and various parameterizations of the set of density matrices were proposed [6]-[12j. see also the 
recent paper [13] for the generalized Bloch vector approach. Recently, two of us proposed another 
parametrization [H] which is based on the particular parametrization of the unitary group found by 
Jarlskog [15] . This parametrization was already applied for modeling quantum gates in the context of 
quantum computing [16j and to investigate possible time evolutions for density matrices [T^. 

In the present paper we propose a similar parametrization for composite systems. This problem 
is of significant importance for analyzing quantum entanglement. Although the analysis of states is 
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parametrization invariant (as much of the problems in physics) a clever parametrization may consid- 
erably simplify the problem and sheds new light into the structure of quantum states. Moreover, as 
usual, the particular parametrization depends very much upon the problem one would like to analyze. 

Recently, a simple parametrization enabling one to analyze the PPT property (positive partial 
transpose) was proposed in [T7|. This special parametrization was used to define a new class of states 
(so called SPPT states) which are PPT but of a very special form. It was conjectured that SPPT 
states are separable. Now, we study the problem of parametrization of composite systems from a more 
general perspective. It is expected that the proposed parametrization will be helpful in analyzing the 
intricate structure of quantum entanglement. 

The paper is organized as follows: in Section[2]we recall the basic ingredients of the parametrization 
of n-level systems used in [T3] . Then in Section [3] it is shown how to generalize it for arbitrary n®m 
composed systems. The procedure is illustrated by explicit examples where it is easy to check for 
separability. Finally, we end with some conclusions. 



2 Parametrization for n-level system 

The parametrization used in [JL4j is defined as follows: any density matrix /j„ may be written as 

Pn = UnD{\i,...,\n)Ul , (1) 

where the matrix -D(Ai, . . . , A„) is diagonal with on the main diagonal and [/„ is a unitary matrix 
from SU{n). Now, following [15j any element Un from SU{n) may be factorized as follows 



where 



(3) 

with Xj E su{n) for j = 1, 2, . . . , n. The antihermitian traceless matrices Xj entering ([3]) are defined 
as follows: Xi is diagonal and the Xj, ioi j = 2, . . . ,n, are given by 





f o.-i 







Xj = 












{ 








(4) 



where denotes k x k null matrix, and \zj) denotes a complex vector from ^ 

( ^ij 



(5) 



together with {zj\ = {zij, . . . , Zj_ij). Taking into account that A\ is diagonal, formula ^ implies 



A-n ■ ■ ■ A^n -^(-^1) • • • ) ^n) A'^ . . . A^ . 



(6) 



Hence, to parameterize p„ one needs (n — 1) complex vectors Z2, ■ ■ ■ ,Zn, with Zj GO ^, and (n — 1) 
real parameters Ai, . . . , A„_i. All together (n^ — 1) real parameters (to parameterize n x n Hermitian 



2 



matrix one needs indeed r? real parameters - n on the diagonal and n{n — 1) off-diagonal - and the 
normalization eliminates one parameter). Now, a simple calculation [16j gives 



V In-j / 

with Ifc being the k x k unit matrix and Vn being the j x j unitary matrix given by 

yj ^ f h-l-(.^-CjWj){zj\ Sj\zj) \ 



where \zj) denotes the unit vector 



that is, 
and 



-Sj{Zji 



(9) 



{zj\zj) = = 1 , (10) 

(11) 



with 

9j := \\zj\\ . 

As usual, in Eq. ([8]) \zj){zj\ denotes the {j — 1) x (j — 1) matrix defined by 



(12) 

i\zj){zj\)ki ■■= ZkZi , (13) 



for A;, / = 1, . . . , j — 1. 

Therefore, pn is parameterized by {n — 1) eigenvalues Ai > ... > A^-i, {n — 1) unit vectors 
i2, • • • ,^n, i-e., Zj defines a point on the {2j — 1) -dimensional unit sphere S'^-'"^, and (n — 1) angles 
02, ■ ■ ■ ,0n from the hyperoctant of (n — 1) -dimensional space. All together we have the correct number 
of independent parameters of a n x n density matrix 

n-l 

J2{2j-l)+2{n-l)=n^-l . 

As an example consider the simplest system, i.e., a qubit corresponding to n = 2. One obtains [13] 

4 . v.- ( f ) . ,14, 

with \z) = e**^ and c = cos 9, s = sin 9. One has therefore 

c^Ai + s^Aa sce^^(Ai-A2) A 
.ce-^Ai - A2) c^X2 + s^Xi J ■ ^ ' 

The above parametrization reproduces the standard Bloch ball: ip and = 29 are nothing but the 
spherical angles on the unit Bloch sphere, and r = Ai — A2 G [0, 1] is the radial coordinate. For Ai = 1 
one obtains the celebrated Bloch sphere of pure states 
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COS 9 
e**^ sin ( 



(17) 



that is 

m = 

up to an overall phase factor. 

3 Composite n<S)m systems 

Consider now a density operator for the composite system hving in ® C™". It is clear that we may 
parameterize it as a density operator living in C"'™'. However, this way we loose information about 
the particular tensor product structure of the total Hilbert space C"™". To control the division into 
subsystems C""* = C" let us consider p as an n x n matrix with m x m blocks, i.e. 



Pn,m = ^ \i){j\(^Pij , 



(18) 



with pij being mxm complex matrices. Our aim is to provide a suitable parametrization for positive 
block matrices. Let -D(Ai, . . . , \nm) denote a diagonal nm x nm matrix with Aj > and Aj = 1. 
It is clear that 

Pn,m = Un,m ' D{Xi, . . . , Xnm) ' Ul^^ , (19) 

where Un^m S SU{nm). Any special unitary matrix Un^m may be written as 

Un,m = , (20) 

where X is an nm x nm anti-hermitian matrix and hence it may be represented as follows 

X = Xi+X2 + ...+Xn , (21) 

where Xi is anti-hermitian block-diagonal and Xj for j > 2 are n x n block anti-hermitian matrices 
with mxm blocks defined as follows: 



X, 





\Zj) 





\ 















la-j ( 





(22) 



where, instead of (n — 1) column vectors zj we take (n — 1) column block vectors 



with Zi j being mxm matrices. Similarly 



(23) 



{Zj\ — {zlj, zlj, ... , ^j_ij) • 

Using the parametrization of block anti-hermitian matrices ()2ip and (I22p we are ready to define the 
following parametrization of the unitary group: 



(24) 
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Hence, we consider unitary matrices from SU (nm) of the following form 



(25) 



where An^m are unitary block matrices (j24p and ^^'j^ is unitary block diagonal, i.e.. 



n,m 



\ 

































Un I 



(26) 



where Uk are m x m unitary matrices. It is clear that 

is a positive block diagonal matrix. Let us denote it by i?(Ai| . . . |A„), where stand for m x m 
diagonal positive blocks 

Afc = UkD{\km, ■ ■ ■ , ^km+m-lWl . (27) 



If m = 1 one has A^ = > 0. Moreover, we add the normalization condition 

Tr(Ai + ... + A„) = 1 . 

One has finally 



Pn,m — ^n,m • • • ^n,m-^(Al| • • • |An)^n|m ■ • • 



t 



(28) 
(29) 



To apply the above formula one needs the explicit form of the unitary components An,m- A straight- 
forward calculation gives (we follow [16]) 



V 



\ 



(30) 



Im j 



where the unit block 1^ appears n — j times. In the above formula Vn,m is a j x j block unitary matrix 
with m X m blocks defined as follows: 



(g)Irrt - \Zj) 


<S) (Irn - Cj) 




\Zj)Sj 




Sj{Zj\ 







(31) 



and \Zj) denotes the normalized block vectors, that is, 

Zk,j 



Z, 



where 



(32) 
(33) 
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Moreover, \Zj 



(Im — Cj) {Zj\ stands for the following (j — 1) x {j — 1) block matrix 



\k){l\®zlf,Zi, , 

kl=l 



and 
with 



Cj = cos , Sj = sin ! 



(34) 

(35) 
(36) 



It is clear that for m = 1 we recover the parametrization used in [14j . 

4 Examples 

Class 1 . Let us consider a 2 (g) 2 system to illustrate our parametrization for the well known 2-qubit 
states. Taking 



Ai = ©2 , A2 = ^(l2 - cr^) , 5' = sinQ;][2, C = cosa][2, U = ax , 
one obtains a family of rank-1 projectors 



P(a) 



(37) 



sin^ a 








sin a cos a 


























sin Q cos a 








cos^ a 



(38) 



which corresponds to a pure state 



tpa = sin a |00) + cos a 1 11) . 



Note, that this state is separable if and only if 5" = or C = 0. For S = C = one ob- 

tains a maximally entangled state. It shows that a nontrivial rotation by a does produce quantum 
entanglement. As a second example in this class let us take S = C = l2/\/2, U = ax and 



Ai 



1 / 1 -p 
4 I 1-p 



A2 



1 / 1 -p 
4 I l + 3p 



(39) 



with — l/3<p<lto guarantee positivity of the matrices A. One obtains the following 1-parameter 
family of 2-qubit states 



m = 4 



/ i+p 
i-p 








2p \ 




(40) 



1-p 
\ 2p l+p J 

This is the well known family of isotropic states which is known to be separable if and only if p < 1/3 
Actually, a point p = 1/3 is not distinguished by our parametrization. 
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Taking S = sin a I2 and C = cos a I2 one obtains a more general 2-parameter family 



I{p, a) = ^—^ I2 ® I2 + pP{a) 



(41) 



which is separable if and only if 



< 



l + 2sin(2a) 

The above example may be generalized as follows. Instead of p9|) let us consider 

P3 



with Pfc > and pi + P2 + Ps + P4 = 1- Taking 



Pi 



(42) 



sin a 
sin/? 

and U = (Tx one obtains the following family 



cos a 
cos/? 



a,/3e [0,7r/2] , 



/ Picl+P2sl 











{P3-PA)SaCa 
V (Pl -P2)Sf3Cf3 



(Pl -p2)Sf3Ci3 \ 

iP3-P4.)SaCa 











(43) 



where 



Srv = sm a , Cn = cos a 



and similarly for S/3,C/3. We stress that one has p{p;a,f3) > and Trp(p;a,/3) = 1 for any a, (3 and 
the arbitrary probability vector p = {pi,P2,P3,Pa) by construction. Interestingly, the above family 
belongs to the class of 2Cg) 2 circulant states considered in [18j. The Peres PPT criterion [19] gives the 
following separability conditions 



P34 + 



> 



\Pl -P2\Sf3Ci3 , 



Pl(?a+P2si > \p3-Pi\SaCa- 

For a = P = Tr/4: the above family reduces to the family of Bell diagonal states 



(44) 
(45) 



( P1+P2 



P(P) 







P3+Pi 



P3-Pi 
\PI-P2 



Pl-P2\ 
P3-PA 



P3+P4. 

P1+P2 J 



(46) 



Moreover, separability conditions (j44p - (j45p reduce to pk < 1/2 for k = 1,2,3,4. Note, that even if 
/o(p) is entangled /?(p; a, /3) might be separable. Consider e.g. Pi = P2 = Pa = 1/8 and p4 = 5/8, that 
is, p(p) is entangled. Now, (I44p is trivially satisfied and (I45p implies sin2Q < 1/2. Hence, /9(p;Q,/?) 
is separable for a < 7r/12 and arbitrary /?. 



C/ass An arbitrary state of a 2 m system corresponds to 



^2,m 



2,m 



zczt zs 
-sz^ c 



(47) 
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with Z = U U{m) and again C = C0SH2, S = sinH2. One finds 



P2,- 



U{CU^AiUC + SA2S)W 


U{SA2C -CU^AiUS) 


{CA2S - SWAiUC)W 


CA2C + SWAiUS 



u 




f CU^AiUC + SA2S 


SA2C -CWAiUS \ 




0,n 


Om 


Im J 


V CA2S - SWAiUC 


CA2C + SWAiUS ) 


V o„ 


Im, 



Note that for S = or C = one obtains a class of block-diagonal matrices 



Ai 






A2 



or 



A2 


Om 




Ai 



(48) 



(49) 



respectively. Being block-diagonal these matrices represent separable 2®m states. It shows that 
quantum entanglement arises only for nontrivial H2 corresponding to C 7^ and S 7^ 0. 
Note, that for 

1. Ai = A2 = A , 

2. [A, [/] = , 

one obtains the following class of 2 m states: 



(50) 



u 


Om 


Om 


Im 



A 


B 




A 





Om 




Im 



with 

A = CAC + SAS , 
Now, if UAU^ = A , then one gets 



B = SAC - CAS 



A 


UB ■ 


(UB)^ 


A 



(51) 

These are block Toeplitz positive matrices and it is well known that they are separable [20]. In this 
way we define huge family of bipartite separable states. Another class of separable states is defined 
by block Hankel positive matrices [20]: taking [/, Ai, A2 and S2 satisfying 



1. [WAiU,E2]=0, 

2. [A2,H2] = , 

one obtains the following class of 2 m states: 



with 



Now, if 



u 


Om 


Om 


Im 




B' 





B' 


B' 


A2 









Im , 



(52) 



UB' = B'U^ 



8 



then one gets 

/ TT A.TT^ Y \ 

(53) 



' UAiU^ 


X 


X 


A2 



with X := UB' . These are block Hankel positive matrices and hence separable [20]. 

Class 3. An interesting class of bipartite n®m states corresponds to Ai = . . . = A„_i = Om and 
An = ^Im- In this case one has 



m 



Pn,m = ^Al^ D{Om\ . . . |0„|I^) All^ . (54) 

This class is parameterized by a positive matrix H„ and (n — 1) complex matrices Zi^n, ■ ■ ■ , Zn-i,n 
satisfying Ylk=i ^1 n^k,n = Im- Taking into account a polar decomposition Z^^n = PkUk, with positive 
Pk and unitary Uk, one replaces the above constraint by 

+ . . . + P^_, = I„ . (55) 

The above equation defines a nonabelian sphere [21j. Therefore, to parameterize the class defined by 
dSlD one may use (n — 1) unitaries {Ui, . . . , Un-i} and (n — 1) positive operators {Pi, . . . , P„_2, H^} 
(since Pn-i niay be calculated from ([55])). Note, that for fixed {Pi, . . . , Pn~2,'^n} one obtains a 
(n — l)m^-dimensional subspace which may be called a nonabelian (n — l)-torus. It is, therefore, clear 
that a class (|54p generalizes the subspace of pure states for n-level single system and hence it may be 
regarded as a nonabelian generalization of the complex projective space CP'"-i [HIS]. Let us observe 
that if Zfc j are normal, that is 

^kjZl- = Zl-Zkj , (56) 

then ()54p defines a rank-m projector which generalizes rank-1 projector (a pure state) for a single 
n-level system. In particular for n = 2 one generalizes a 2-dimensional Bloch sphere (one may call it 
a nonabelian Bloch sphere): 

1 



52 


use 


csw 





(57) 



which is parameterized by two nonabelian angles: unitary U and positive H2 (C = C0SH2, S = sinS2). 
All together 2m^ parameters. Note that for n = m the class (|54p defines a set of extremal states of 
the extended quantum theory proposed recently by Zyczkowski |22j . 



5 Conclusions 

We proposed a parametrization of density matrices of composed n®m quantum system. For m = 1 
this parametrization reduces to the one used recently in [Ijj. Note, that it may be generalized for 
multipartite systems living in C"^ ® . . . ® C""^ . Indeed, instead of dealing with n x n block matrices 
with m X m blocks, in the multipartite case one has to consider rii x ni block matrices with blocks 
being 722 x 722 block matrices with blocks being 713 x block matrices and so on. Although the 
strategy seems to be simple the technical part of the story is quite involved. It is anticipated that 
the presented parametrization will find interesting applications in the study of quantum properties of 
many partite systems. 
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